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Abstract 

A modified Borel transform of the Adler function is used to resum the 
hadronic tau decay width ratio. In contrast with the ordinary Borel trans- 
form, the integrand of the Borel integral is renormalization-scale invariant. 
We use an ansatz which explicitly accounts for the structure of the leading 
infrared renormalon. Further, we use judiciously chosen conformal transfor- 
mations for the Borel variable, in order to map sufficiently away from the 
origin the other ultraviolet and infrared renormalon singularities. In ad- 
dition, we apply Pade approximants for the corresponding truncated per- 
turbation series of the modified Borel transform, in order to further ac- 
celerate the convergence. Comparing the results with the presently avail- 
able experimental data on the tau hadronic decay width ratio, we obtain 
a s (M 2 ) = 0. 1 192 ±0.0007 exp .±0.0010 E w+CKM± 0.0009th. ±0.0003 evo i.. These 
predictions virtually agree with those of our previous resummations where we 
used ordinary Borel transforms instead. 
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I. INTRODUCTION 



Extensive perturbative calculations in QCD have made available the truncated pertur- 
bation series (TPS) of various observables to the next-to-next-to-leading order (NNLO), 
i.e., including the terms ~a^. There is a longstanding problem of how to extract (resum), 
in a reasonable manner, the values of such observables S(a s ) as functions of the strong 
coupling parameter (M z 2 ). Many of the resummation methods are based solely on the 
available TPS SWai^). Some of these methods eliminate the (unphysical) renormalizat ion- 
scale (RScl) and renormalization-scheme (RSch) dependence of the TPS by fixing the RScl 
and RSch in the TPS itself. Among them are: the RScl fixing of Brodsky, Lepage and 
Mackenzie (BLM), motivated by large-n/ arguments |T|; the principle of minimal sensitivity 
(PMS) 0; and effective charge method (ECH) Some of the more recent methods are: 

the approaches using "commensurate scale relations" |J; ECH-related approaches [fTHlCfl; a 



method using an analytic form of the coupling parameter ||11|| ; a method using expansions in 
the two-loop coupling parameter |f[2]j ; a method which disentangles the running coupling and 
conformal effects (skeleton coefficients) [O] ; methods using conformal transformations either 



for the coupling parameter |14|] or for the Borel expansion parameter [T^|T^]; the method 
of Pade approximants (PA's) [|T7|-|TP||; an RScl-invariant [^DJ and RScl- and RSch-invariant 
extension of the PA-approach [[21|j . 

Some of the aforementioned approaches allow for incorporation of additional, nonpertur- 
bative, information on the observable, e.g., the information on the location of the dominant 
infrared (IR) renormalon singularity pl| via the fixing of the remaining free parameter in 
the approximant, or the information on the location of various renormalon singularities via 
leading-/^ resummation |||| . Full information on the known leading IR renormalon struc- 
ture can be incorporated in a natural way through an explicit ansatz, in the approaches 
using Borel transformations of the considered or of the associated observable (see, e.g., 
|2"2|). Incorporation of the available leading renormalon information appears to influence 
significantly the numerical predictions, particularly for low-energy (low-Q 2 ) QCD observ- 
ables, through a more accurate description of the Borel amplitude in the region between 
the origin and the leading IR renormalon singularity. This method has an advantage over 
the widely used (unresummed) Operator Product Expansion (OPE) approach [p3 |, in which 



the perturbative contributions are usually taken as TPS, and additional nonperturbative 
terms oc 1/Q 2n appear (cf. [^4j and references therein). Since the perturbation series is 
divergent, the latter approach has the problem that there is no natural separation between 
perturbative contributions and power-suppressed nonperturbative terms, a problem which 
does not occur in the former method. However, the Borel-resummed value for an observ- 
able, with the renormalon information incorporated in it (e.g., via an explicit ansatz), in 
principle still does not represent the full value of the observable, because power-suppressed 
terms should be added to it. Dispersive models ( |^5| and references therein) contain and 
predict the OPE-type power-suppressed terms. Recently, it has been argued |26| that the 
power-suppressed terms can be obtained from the knowledge of the perturbation series, or 
parts thereof, and of the IR renormalon structure of the observable. 

In all these resummation approaches, for obvious reasons, it is highly preferable to con- 
sider QCD observables whose values are presently known at a high precision level. Presently, 
high precision experimental data on r lepton decay widths p7|-p2| are available. The non- 
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strange hadronic r decay width ratio R T (AS = 0) can be obtained by using constrained 
fit values |33j of the (basis modes) leptonic branching ratios B e = B(r~ — > e~V e v T ) = 
(17.83 ± 0.06) x 10~ 2 and = B(t~ -> yTV^Vr) = (17.37 ± 0.07) x 10~ 2 of r, and sub- 
tracting the strangeness-changing part |32j R T (AS 7^0) = 0.1630 ± 0.0057 

flr( As = o) B r( i: r ^ M i°7y » - MAsm (i) 

r(r -> v T e v e {l)) 



(l-B e - Bjj) 
B P 



R T (AS^0) = 3.4713 ±0.0171 . (2) 



The above ratio is a QCD observable at relatively low momenta y/q 1 ~ m T w 1.8 GeV, 
experimentally known to high precision, thus presenting an experimental challenge to the 
theory. The challenge consists in predicting the strong coupling constant a^ s (M|) so that 
the theoretical uncertainty (Sa s ) t h.., which partly originates from the uncertainty of the 
method of resummation and partly from the uncertainty of the associated Adler function, is 
smaller or comparable to the uncertainty (Sa s ) exp , originating from the small experimental 
uncertainty SR T ~ ±0.017 given above. 

There is now a wealth of theoretical results [p4| , p4| -|40|| available on the observable ([|), 
and, in particular, on the associated Adler function D(Q 2 ) - perturbative as well as nonper- 
turbative. The main QCD theoretical problem connected with the observable (|2|) is that the 
pertaining momenta are somewhat low \Q\ ~m T (~ 1.8 GeV), and so the relevant coupling 
parameter a s (Q 2 ) in the perturbation expansion is large. It is thus important to take into 
account, in any resummation procedure for R T and/or D(Q 2 ), not just the known pertur- 
bative coefficients, but also a significant part of the nonperturbative information, i.e., the 
leading infrared renormalon. On the other hand, the nonperturbative contributions to R T 
that are represented by power-suppressed OPE terms (apart from the well known quark 
mass contributions) have been shown to be consistent with zero in the ALEPH analysis J27| . 



One version of this program was carried out in our previous work p2| . There we 
first used the known information on the leading infrared (IR) renormalon and on the 
corresponding l/Q 2 -suppressed term of the OPE for the Adler function D(Q 2 ), in or- 
der to predict the 0(aj) coefficient d?,{Q 2 ) of D{Q 2 ). Using this prediction, and ju- 
diciously chosen conformal transformations, we resummed R T by employing an ordinary 
Borel transform D(b) of the Adler function with an ansatz which explicitly accounted for 
the leading IR renormalon structure. Comparing the obtained expression with the ex- 
perimental values (0), we obtained definite predictions for the strong coupling parameter 
a s (M 2 ) = 0.1193 ± 0.0007 cxp . ± 0.0010 ew +ckm ± 0.0009 mcth . ± 0.0003 cvoL . We fixed the RScl 
by the principle of minimal sensitivity (PMS) applied to the approximant, and by choosing 
the MS RSch. 

In the present work instead, in order to obtain a cross-check of the predictions, we em- 
ploy not ordinary Borel transforms but modified Borel transforms of D(Q 2 ). Such Borel 



transforms were introduced by Grunberg [fil|j , on the basis of the modified Borel transforms 
of Ref. [ 4"2]| . One attractive feature of Grunberg's transforms is that the integrand in the 
Borel integral is RScl-independent, in contrast to the case of the ordinary Borel transforms. 
Another interesting feature is that they represent integral transformations of a significantly 
different form than the ordinary Borel transformation, and have therefore a different singu- 
larity structure. Therefore, their application to the hadronic tau decay width ratio and the 
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subsequent extraction of the prediction for a s (M 2 ) could represent a powerful cross-check 
of the results of previous work |22[] based on ordinary Borel transformations. 

In Sec. [Tl] we recall the known basic theoretical formulas for R T (AS = 0) and the associ- 
ated Adler functions, as well as the reduction to the massless QCD observable r T . In Sec. |T|. 
we present the modified Borel transforms D(b) of the massless Adler function D(Q 2 ). We 
refer to the Appendix for details about D(b). In Sec. |IV| we perform the resummation for r T , 
by contour integration of D(Q 2 ) in the complex momentum plane, explicitly accounting for 
the leading IR renormalon structure of D(b) through an ansatz, choosing conformal trans- 
formations of the Borel variable b to map away other singularities of D(b), and employing 
Pade approximants. In Sec. [V] we compare the obtained expressions with the experimental 
results and determine a s (m 2 ) and a s (M 2 ). We further estimate the theoretical uncertainties 
of the prediction. Sec. |VT] contains a summary and brief discussion of the differences between 
our results and those of other analyses of R T . 



II. THE KNOWN BASIC FORMULAS, REDUCTION TO MASSLESS QCD 



The restriction AS = in Eq. means that only hadrons with quarks u and d are 
produced. Thus the observable (H) is already close to being massless. This fact removes 
some of the complications in the theoretical analysis. 

The ratio (0) can be expressed, via the application of a variant of the optical theorem, 
and the subsequent use of Cauchy's theorem and integration by parts, as a contour integral 
in the complex momentum plane (see, for example, Refs. piMOj): 



{AS = 0) = 



R T {AS = 0) 
3\V ud \ 2 (l + 5 EW ) 



ew; 



-7T1 



7 1 



3 r 



mi 



1 + — )D 

mi 



LfT/ 



-s) + -D L (-s) 



- 1 . 



(3) 
(4) 



Here we factored out, for convenience, the square of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element \V u d\, the electroweak (EW) correction parameter 5 EW = 0.0194 ± 
0.0050 ||I3|, 4"4] , |24 |, and the residual EW correction parameter <5^ w = 0.0010 ^3|. The con- 



tour integration in Eq. (Q) is counterclockwise in the complex s-plane, and the general Adler 
functions D 1 ^ and D h are related to the V+A current-current correlation functions as: 



D 



LfT, 



d x 

J=0,1 



ds 
mi ds 



n 



(0) 



(0) 
ud,V 



n 



ud.A\ 



Here, J is the spin of the hadronic system in its rest frame (L: J = 0; T: J=l); ^ U dv/A 
the components in the Lorentz decomposition 



(5) 
(6) 
are 



of the two-point correlation functions ^dv/A °f ^ ne vec t° r K!d 
Aud = dl^l5U (color-singlet) currents 



L ud,V/A 



ad V/ AVI ! (7) 
d^u and axial-vector 
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- in 



ud,v(l) 



Ui ud,A 



in). /^e^-(0|T{< d (x)^(0)t}|0) 



(8) 
(9) 



In the massless quark limit {m u ,d ~^ 0), D L (s) vanishes and the vector and axial-vector 
contributions to D LrfT become equal in perturbation and D LfT (— s) —>■ (1 + D(—s))/(2ir 2 ), 
where D(Q 2 ) is the canonically normalized massless Adler function with the perturbative 
expansion^] 



oo 

1 + d na n 

n=l 



(10) 



Here, a = a s (fi 2 ; c 2 , C3, . . .)/tt is the QCD couplant at the renormalization scale (RScl) /i 2 
and in the renormalization scheme (RSch) characterized by the coefficients Cj (j > 2) in the 
beta function 



^4^ a = P( a ) 



dfi 2 



-(3 a 2 [l + c x a + c 2 a 2 + ■ 



(11) 



Here, (3q = (11 — 2n//3)/4 and c\ = (102 — 38n//3)/(16/?o) are two universal constants which 
depend only on the number of active quark flavors nf. The Adler functions are quasiobserv- 
ables, in the sense that they are independent of the RScl and RSch. 

In order to apply the massless QCD analysis to r T (A5' = 0) (|j), we have to subtract from 
it the quark mass (m u ^ 7^ 0) contributions. This can be carried out p4] within an operator 
product expansion (OPE). The largest quark mass contributions are quark condensate terms 
of dimension d = 4 (cxl/m^) 



5r T (AS = 0) mud -£ ~ 16tt 

167T 2 / 2 m 2 



2 ym u + m d )(qq) 



mi 



1 , 23 ( a s {m 2 T ) 



7T 



mi 



1 , 23 f a s (m 2 T ] 



TV 



-0.00265 x (1 + 0.03) w -0.0027 



(12) 
(13) 



In Eq. (|T2]) we denoted (qq) = (mm) w (dd). Here renormalization scale can be taken to be 
/i ps m r . In Eq. (|13|) we used the PCAC relation (m u + m d )(qq) ps —j 2 m\ (^ = 92.4 ± 0.3 
MeV; m„— = 139.6 MeV). There are corrections to this relation and to expression ( fEf ) of the 
order ~m 2 d /m 2 , i.e., of the order of the OPE d—2 terms which can reach, at most, values 
~ 10~ 4 . In obtaining the numerical value in Eq. flilf ), we further used m T = 1777 MeV and 



a Amt 



MS) ~ 0.32. 



The OPE approach of [2^| includes other power-suppressed nonperturbative terms that 
contribute to r T but do not stem from quark masses: the d = 4 gluon condensate term, and 
the d = 6 term. The latter term could be large, but it has also comparably large uncertainties 



1 The (ud) Adler functions D L+T and D L (S)— (0) usually include, by convention, the additional 
CKM factor \V ud \ 2 (e.g., see Refs. |0|). 
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|24| . The gluon condensate contribution to r T in the OPE is a 2 -suppressed. The ALEPH 



analysis |27j indicates that these d = 4, 6 nonperturbative contributions to r T are consistent 
with the value zero, £r r (NP; m Ut d = 0) = 0.000 ± 0.004. We should keep in mind, however, 
that the ALEPH analysis assumed that the part of the associated Adler function which has 
no power-suppressed terms is a (N 3 LO) TPS, while we will perform resummations of this 
part by taking into account its renormalon singularity structure. Nonetheless, we consider 
the ALEPH analysis as indicative that even in our resummation framework the power- 
suppressed OPE-type terms in ry (apart from the quark mass terms) are either consistent 
with zero or very small. Therefore, we will ignore in our analysis of ry any OPE power- 
suppressed nonperturbative terms other than those in Eqs. (|12|)-(|T3|). 

We thus regard, in our framework, as nonperturbative massless contributions only those 
contributions which appear as a consequence of explicit IR renormalon structure of the 
Borel transforms in the resummation. For example, the leading IR renormalon of the Adler 
function D(Q 2 ) (which we will account for) gives contributions to D(Q 2 ) which can be 
partially represented as a d = 4 power-suppressed term oc l/Q 4 - This, however, does not 
necessarily mean that there is no additional, genuine OPE-type d — A term (oc (aGG) /Q A ) 
in D(Q 2 ). The uncertainties of the OPE d — A, 6 massless terms as given by ALEPH are 
large enough to accommodate the possibility of significant nonzero values of these terms. 
For example, (aGG) = 0.001 ±0.015 GeV 4 . Toward the end of Sec. |V| we will briefly discuss 
how the latter uncertainties would influence the final prediction for a s in our analysis. 

The ALEPH analysis further assumed that massless d = 2 terms (oc 1/Q 2 ) are not present 
in the Adler function (and thus in ry), as suggested by the OPE. Such terms were suggested 
by the authors of Refs. P5lff6j as an effective tachyonic gluon mass contribution reflecting 



nonperturbative short-distance QCD. However, the authors of Ref. [47] showed that the 
strength of such terms is consistent with zero.Q 

When subtracting the quark mass contributions fli~2|)-(ll3D from Eq. @, we end up with 
the massless QCD observable 

ry = r r (A5 = 0;m M = 0) = r T (A5 = 0) - 5r T (AS = 0) mu ^ (14) 

3 / \ 

-s) , (15) 

= — r dy (l + e iy ) 3 (l - e ly ) D(-s = m 2 T e iy ) , (16) 

with the canonically normalized massless Adler function D(Q 2 = —s) defined in Eq. (|TCJ|). 

The perturbation coefficients dj = dj(fi 2 /Q 2 ; c 2 , . . . , Cj) in Eq. ([TOD depend on the RScl 
fi 2 and the RSch parameters Cj in a known specific manner, because D(Q 2 ) is RScl- and 
RSch-independent. Knowing them at a specific RScl and RSch, we know them at any RScl 
and RSch - see Eqs. (|A.5|) - ([A.8 ) in the Appendix. The first two coefficients have already 




been calculated, and in the MS RSch and at RScl jjL 2 = Q 2 , for n/ = 3, they are: df' = 1.6398 



2 They did this by fitting a d = 2 finite energy sum rule, which is apparently well satisfied at the 
used relevant scales [48,49], to the new ALEPH data on spectral functions extracted from the r 
decay measurements. 
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SOfl ; c?2 = 6.3710 [|5lJ (the superscript (0) denotes the value at the aforementioned RScl 



and RSch). 

In order to extract the experimental value of r T ([14]), according to Eq. (|3[) we need to 
use the values R T (AS = 0) of Eq. (|), 5 EW = 0.0194 ± 0.0050 P,g,@], 5' EW = 0.0010 |3 



5r T (AS 0) mu d ^o QT3|), and in addition the value of the CKM matrix element \V U d\, which 
we take as0 

\V ud \ = 0.9749 ±0.0021 . (17) 

This leads us, via Eqs. @, fl3|) and (|13|), to the following values for the massless QCD 
observable (|T4|) 

r T = r T (AS = 0;m Uid = 0) 

= 0.1960 ± 0.0059 exp . ± 0.0059 EW ± 0.0051 C km (18) 
= 0.1960 ±0.0098 . (19) 



In Eq. flT9|), the three uncertainties of Eq. ( P~8| ) were added in quadrature. 

The values (p~8|) , together with the contour integral expression ([H]), will be the starting 
point for our massless QCD resummation analyses of hadronic r decays. 



III. MODIFIED BOREL TRANSFORMS 

The most straightforward way to perform the resummation for r T given in Eqs. ([15])- 
(|T6|) would be to insert the known truncated perturbation series (TPS) for D(—s) given in 
Eq. (|I]) and perform the momentum-contour integration, i.e., the approach of ||39|| . The 
method is, firstly, fraught with ambiguities from the choice of RScl and RSch. The final 
result contains residual, but significant, RScl and RSch dependence, due to the truncation 
of the series in Eq. (|i~0[). Secondly, the method does not incorporate the known renormalon 
structure of the Adler function's Borel transform D(b). 

In a previous paper f22fl , two of us addressed the second problem, by employing in the 
resummation a (ordinary) Borel transform D(b) that includes the leading IR renormalon via 
the ansatz D(b) = R{b)/ (1 — 6/2) 1+I/ , and by introducing in addition conformal transforma- 
tions b = b(w) in order to map sufficiently far away from the origin the singularities of the 
UV (and the remaining IR) renormalons. However, the integrand in the (ordinary) Borel 
integral is RScl- and RSch-dependent. This can be inferred from the definition 



3 The standard model (SM) unitarity-constrained fit predicts \V uc i\ = 0.9749±0.0008 [33]. However, 
the values extracted from the decays of mirror nuclei are lower: \V u d\ = 0.9740 ± 0.0010. This 
extraction has significant theoretical uncertainties (see for further References). The values 
extracted from neutron decays are even lower \ V u d\ = 0.9728±0.0012 ( |33[] and References therein), 
but appear to have smaller theoretical uncertainties. In view of all these considerations, we take 
in our analysis the value range as given in Eq. ([17]). Here, the central value is from the unitarity- 
constrained fit, and the uncertainty is increased so that ( |i~7| ) covers all the values from the decays 
of mirror nuclei and the upper half of the interval from neutron decays. 
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D(Q 2 



1 f°° 



Po Jo 



dbexp 



(3 a(p 2 ; c 2 , . . .) 



D{b-^/Q\c 2 



(20) 



and from the expansion of the transform around the origin 



D(b; p 2 /Q 2 , c 2 , . . .) = 1 + *™ ( » U f d Am^l • • • » O (b_ (21) 



1! 



V/3o/ ^2 n! \A>/ 

For example, direct application of the derivative with respect to the RScl /i 2 gives 



d 



P d\np 2 



"V(ft«)£)(M 



-6/(/3 a), 



/3(a) 



+ 1 



+ 



d 2 



'/3(a) 



+ 1 + -(2c x di + c 2 ) 



/?oa 2 / 3 





(22) 



using the notations of Eq. (|Tl|). Once going beyond the one-loop approximation of the 
RGE-evolution (|TTD, the integrand is RScl-dependent. If we knew the exact expression of 
the integrand, the total integral ( ^0|) would be RScl-independent. However, since we have 
a TPS available for D(b), we are forced to use a TPS for R(b) = (1 — b/2) 1+u D(b). This 
truncation then results in the residual RScl and RSch dependence of the resummed result 
for D{Q 2 ) and thus for r T . In we fixed the RScl-parameter £ 2 = [i 2 /Q 2 according to the 
principle of minimal sensitivity (PMS), i.e., dr T /dC, 2 = 0. The RSch was chosen to be MS. 

In the present paper we apply, instead of the ordinary Borel transform, several variants of 
the modified Borel transform D(b) of the Adler function D(Q 2 ). The integrand will be RScl- 
independent. This Borel transform was introduced for QCD and QED (quasi)observables 
by Grunberg |4l]], who in turn constructed them on the basis of the modified Borel trans- 



formations of Ref. j42H . The integral transformation for D(b) is written in the form 



Po Jo 



dbexp 



{pi{Q 2 )+c)b 

A) 



D(b;c) 



(23) 



Here, -D(6;c) = exp(c b/j3 )D(b; 0) and has no Q 2 -dependence; c is a specific arbitrary 
constant; pi is the first RScl and RSch invariant || of the Adler function 



Pi(Q 2 



-cM/iVQ 2 ) + A) In 



Q 2 

Po In =j , 
A 



(24) 



where A is the universal scale appearing in the Stevenson's equation ||, and A is a scale which 
depends on the observable but is RScl- and RSch-independent. We note that pi(Q 2 ) = 
[l/a^- loop \Q 2 ) +c], where c is a constant. Therefore, D(b; c) reduces to the ordinary Borel 
transform D(b; p 2 /Q 2 = 1, . . .) times the factor exp[(c+c)6//3o], when higher than one-loop 
effects are ignored (large-/?o approximation). 

Since D{Q 2 ) and p\{Q 2 ) in Eq. ( p3|) are RScl- and RSch-independent, so D(b; c) possesses 
no explicit RScl- and RSch-dependence. However, as shown in the Appendix, the modified 
Borel transformation (^) is one in a large class of Borel transformations, where each of them 
corresponds to a specific choice of the RScl and RSch. The transformation (|23|) is obtained 
from this class by choosing the so called modified 't Hooft RSch (mtH RSch: Cj = c{, 



S 



j = 2, 3, . . .). We also show that in this case the constant c is in fact the RScl parameter c = 
/3ohi(£ 2 ), where £ 2 = /i 2 /Q 2 . This special Borel transformation has the remarkable property 
that the integrand in the Borel integral is RScl-independent (c-independent). Therefore, 
we can call it RScl-independent Borel transformation, since the change of the RScl only 
changes the convention of separating the (RScl-invariant) integrand in Eq. ( |23"D into two 
factors. 

Yet another useful property of the modified Borel transform (23) is the following: due 



to the very simple Q 2 -dependence of the integrand in Eq. fl2~3|) , the contour integration (p~6[ ) 
in the complex momentum plane can be performed analytically, leading to a rather simple 
expression for the observable r T , as will be shown in the next section. This is not possible 
if the mtH RSch is abandoned. 

The coefficients of expansion of D(b; c) in Eq. ( p3|) around 6 = are related with those of 
the expansion of D(Q 2 ) most easily when we use for the latter expansion the specific mtH 
RSch c fc = c\ (k > 2) and RScl /i 2 = Q 2 



^ M/A)) Clb/ftl f{c~di)b\ \ A , (Jx-d) 
D{b;c) = — exp 




r(l + Cl 6/A)) "V A) [ (l+cxb/Po) \Po 



+ ^ 2 (l + c 1 b/f3 )(2+c 1 b/(3 )---(n+c 1 b/f3 ) Go) J' (25) 

Here, <i,- = dj([i 2 /Q 2 ; c 2 , . . . , Cj) with /x 2 = Q 2 , = (fc = 2, We refer to the 

Appendix for details. Expanding each term on the right-hand side in powers of b then 
yields the expansion of D(b\ c) around 6 = 0. Here we can see again that in the large-/?o 



approximation (ci — > +0) the above expansion reduces to the expansion (21) of the ordinary 
Borel transform D(b; 1, . . .) times the factor exp[(c— di)b/ (3q\. 

The ordinary Borel transform -D(6) is known to have singularities at b= — 1, —2, . . . (UV 
renormalons) and at 6 = 2, 3, . . . (IR renormalons) . The renormalon resummation of D(Q 2 ) 
and of the hadronic r decay width in the large-/?o limit has been performed in Refs. |52]-|54|. 
The IR renormalons on the contour of the Borel integration ( |20|) cause ambiguities, above 
all the leading IR renormalon (6 = 2). The singularity of D(b) at 6 = n has the form 
1/(1 -6/n) 2+7 " +nCl / /3 ° for n > 3, and 1/(1 -b/2) 1+ ^ +2ci / /3 ° at 6 = 2, where i n is the one-loop 
anomalous dimension of the operators corresponding to the l/Q 2n -terms in the operator 
product expansion (OPE) of the Adler function D(Q 2 ). In the case n — 2, it is known 
that y = |]55|] . Thus, at the leading IR renormalon we have .0(6) oc 1/(1 — b/2) 1+u , with 
!/ = 2ci//5 (=1.580, when n/ = 3). 

On the other hand, the modified Borel transform D{b) in Eq. ( J23J has the (IR and UV) 
renormalon singularities at the same locations as D(b), but with simpler powers [f4T| , [42"| 
the IR renormalon singularities are of the form 1/(1 - b/n) 2+ < for n > 3, and 1/(1 - 6/2) 1 
for the leading IR singularity^ [see Eq. ( |A.13| ) in the Appendix, with K = n]. Therefore, we 
can define a new function R(b; c) 



4 This is true even when the mtH RSch is abandoned and the integral transformation (cf. Ap- 
pendix) becomes considerably more complicated - as follows from the considerations of Ref. E3] . 
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R(b;c) = (l-b/2)D(b;c) , (26) 

which has a considerably softened singularity at 6 = 2 (cut instead of pole). This function 
can be determined by resummations of the TPS of R(b;c) (cf. also Ref. fL5f ). The latter 
TPS is known to the same order as D(b; c) - see Eq. (H): 4 0) and (<S> d u d 2 ) are known 
exactly, and djj J3) is known approximately. The same is true for the modified functions 
Rjfac) (j = l, 2, 3)_ discussed below, Eqs. 

The function D(b; c) as defined in Eq. (p3|) has additional singularities, as seen from the 
expansion fl25|) : 1.) it is nonanalytic at 6 = 0, due to the (finite) factor (ci6//3 ) Clb///3 °; 2.) it 
has additional (spurious) poles at b = —(3 /ci (~ —1.266 when n/ = 3), b = —2/3 /ci, etc. In 
approximate numerical evaluations of D(Q 2 ) and r r via -D(6; c), such singularities may have 
disturbing effects unless they are sufficiently far away from the origin. Therefore, we can 
define the following variants which will be used in (approximate) resummations 

R^b; c) = (1 - 6/2)( Cl 6//3 r cib//3( T(l + c 1 6//? )e Clfe//3 °D(6; c) (27) 
= (l-V2)exp[(c + c 1 - rf >//3.]{l+ (A) + . . . } , (28) 

ifefc c) = ^(6; c) (3 + Ci6//5o)2 • (30) 

We recall that Rj(b;c) = exp(c b/fto)Rj(b; 0). In the above expressions, we restored an- 
alyticity at b = by factoring out, instead of the factor (ci6//9o) Clfe// ^°, the combination 
(c 1 6//3 ) Clb//3o e- Clfc / /3 V r ( 1 + cib/(3 ). Our main motivation for this lies in the following: the 
factor (ci&//3o) Clb// ^° is increasing extremely fast with increasing b, while the aforementioned 
combination has weak 6-dependence 

( Cl 6//9o) clWo e- cl ^° 1 . , , Q1 , 

when b — > 00 . (31) 



r(i + ah/fa) ^nab/fa 

Hence, functions Rj(b; c) behave at large b roughly as Rib] c), or D(b; c), 



Rj(b; c) ~ (1 - 6/2)^/27^1 6/ A) £(&; c) when 6 -> 00 , (32) 

i.e., they neither decrease nor increase violently. Therefore, any approximate resummation 
method will have better chance when applied to them than to an extremely fast increasing 
or decreasing version. On the other hand, if we just factored out the factor (ci&/ / 3o) Cl ^ /3 °, 
the resulting function, though analytic at b = 0, would at large b decrease violently, as 
~( Cl b/f3 )- c ^D(b;c). 

The function Ri(b;c) has spurious (unphysical) poles at b = — /3 /ci, — 2/3 /ci, . . . (~ 
—1.261, —2.53, when rif = 3). But the function R2 has a possible advantage over R± in 
resummations, since it has no spurious (nonphysical) pole at b=— (3q/ci, and R3 has no such 
poles at b = —(3 /ci and — 2/3 /ci. 

Since the coefficients and d^ (=> d\, dz) of the (massless) Adler function are known, 
the power expansion of Rj(b;c) is known to the next-to-next-to-leading order (NNLO, in- 



cluding the term ~6 ) via Eqs. (|25|)-(p0|). In a previous work p2 |, we presented an argument, 
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via a bilocal expansion of the Borel amplitude D(b), that d 3 ~ 25 ± 5. We also discussed 
there the estimates of d% presented by other authors, and concluded that the following 
estimate is rather safe 

25 ±10. (33) 

We will use the above values which allow us to obtain the power series of Rj(b;c) up to 
N 3 LO (~o 3 ). 



d 



(0) 



d 3 (p 2 = Q 2 ;MS) 



IV. RESUMMATION PROCEDURE 



We will apply summations to the N 3 LO truncated power series (TPS) of the functions 
Rj(b; c). However, in order to obtain the (resummed) values of the massless QCD observable 
r T (0), we have to perform first the complex momentum contour integration (fl6|) , with the 
massless Adler function there having the integral form (^3|) in terms of the invariant Borel 
transform D(b; Q 2 = m 2 e iy ), i.e., in terms of the related functions Rj (f2~7|)-([30l). This angular 
^-integration can be performed exactly, because the ^/-dependence of the invariant pi (|24| 
appearing in Eq. ( ^3"D is simple 



pi(m^e 12/ ) 



pi(m 2 T ) 



(34) 




{p x {m 2 T )+c) b 



sin(7r&) 




6(1 — 6/3)(l — 6/4) 



D(b;c) 



(35) 



sin(7ro) 



6(l-6)(l-6/2)(l-6/3)(l-6/4) 



(ci6//?o) Cl6//3o e- clf>//3 ° 
T(l + Cl b/p ) 



ifj = l,' 

fj{b) = { (2+c 1 b/p )/(l+c 1 b/p Q ) if j = 2 , 

(3+ci6/A)) 2 /(l+ci6/A))/(2+ci6/A)) if j = 3 . 



where 



(36) 



(37) 



The integration contour in Eqs. (|35|) (p6|) is chosen slightly below (or above) the positive 
axis, in order to avoid possible singularities of the integrand on the positive axis. When 



2. 



knowing reasonably well the RScl- and RSch-invariant functions Rj(b) in Eq. (|3q) , e.g. 
via resummations as discussed below, the massless QCD observable r T becomes an expres- 
sion whose value depends uniquely on the value of the QCD coupling parameter, e.g. on 
MS)=7rao- This dependence originates from the known dependence of the invariant 



a Am. 



pi(m 2 ) in the exponent in Eq. (|36[), the latter being determined by an integrated version of 
the RGE, called also the (unsubtracted) Stevenson equation 



pi(m T ;a J 



4 0) +A)ln 



ml 



A 



4° ) + - + c 1 ln(-^) 
ao \1+Cia y 



dx 



1 



x 2 (l+dx) /^(x 



(38) 
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where the perturbation expansion of the last term is given via the definition ([TTJ) of the 
/^-function 



fe(x)//3 = -x i (l+c 1 x+cfV+c!fV + - • •) . (39) 

Here, the MS coefficients are functions of the number of active quark flavors rif and are 
known up to the N 3 LO (df 8 , four-loop) f56|. The number of active quark flavors is assumed 



to be here nj = 3, because the scale of the process is y\Q 2 \ =m T (^1.777 GeV). 

We can note from Eq. (|36]) that ambiguity of integration over 6 in r T at the first infrared 
renormalon singularity of D(b) at 6 = 2 is suppressed, because the factor sin(7r6) has a zero 
there. This is similar to the one-loop approximation in the approach with the ordinary 



Borel transform |Tj||57|| . However, here the absence of the ambiguity is not due to an 
approximation, it is exact and due to the discussed RScl-invariant Borel approach. 

Although the ambiguity at 6 = 2 is suppressed in the Borel-type of integration fl36|) , we 
wish to emphasize that it is nonetheless very important to factor out the leading renormalon 
factor 1/(1 — 6/2) there [according to Eq. fl26|)1. This is so because any resummation of a 
TPS (6) represents also a quasianalytic continuation of the corresponding function into the 
region away from the origin, and such continuation is of a better quality when there are as 
few singularities near the origin as possible. The functions which we will resum in Eq. ([36]) 
are Rj(b;c), i.e., the functions which have their original pole singularity at 6 = 2 factored 
away according to Eq. (|26|). The remaining singularity in Rj(b; c) is then significantly weaker, 
it is a cut of the type ~ ln(l — 6/2) [see Eq. ( A.13|) with a = 1 and k = 2, in conjunction 



with Eq. (|A.20|) ] . In this way, in the relatively wide region 6 ~ 2 we can achieve reasonably 
good values of the integrand of Eq. (|36"D, which leads then to good predictions for r T via 
Eq. ([3(J) as functions of a^ s (m^). If we did not factor out the mentioned singularity, the 
obtained (resummed) values would definitely lead to bad values of the integrand at 6~2, the 
fact which would influence the predicted values of r T and thus of a^ s (m^). The numerical 
importance of factoring out the leading IR renormalon singularity of the Borel transforms 
of the Adler function at 6 = 2 before doing resummation was pointed out in our previous 



work |22j| , which involved ordinary Borel transforms. There we showed that the predicted 
values of a^ s (m^) from ry depend crucially on whether this factoring-out procedure has 
been performed. 

The functions Rj(b; c) in Eq. fl36]), whose TPS's we want to resum via methods of quasi- 
analytic continuation, have the singularities near the origin which are negative (at 6 = — 1 
and lower). The nearest singularity on the positive axis is at 6 = 3. At 6 = 2 (IR2) there 
is a weaker logarithmic singularity. The negative singularities near the origin constrain the 
convergence radius of the perturbation (power) series of Rj(b]d)'s to r = l, and thus repre- 



sent a possible hindering element to efficient resummations. In our previous work [22|] we 
proposed how to extend the convergence radius up to IR2 - by either of the following two 
conformal transformations w = w(b): 



VTTb- Jl-b/3 VTTb - Jl - 6/4 

w 3 (b) = V / = , w A {b) = V , = . (40) 

VTTb+Jl- 6/3 VTTb + Jl- 6/4 



Conformal transformation 1^3(6) maps all the renormalon singularities to the unit circle in the 
w-plane, except for the first IR renormalon which is mapped to ^3(6 = 2) = 1/2. We further 
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note that iu 3 (& = 3) = 1, and u> 3 (5= — 1) = — 1. Furthermore, all the spurious singularities 
b = —n(3 /ci (< —1.265) are also sent to the unit circle. Conformal transformation w^(b) 
does the same thing, except that the first two IR renormalon singularities are both inside 
the unit circle: w 4 (6 = 2) ps 0.42, u> 4 (5 = 3) = 0.6. We further note that w i {b = A) = 1, and 
Wi(b — — 1) = —1. The inverse transformations are 

w s 3w 3 (16/5)w 4 , . 

(l-w 3 + w£) (1 - (6/5)w 4 + W4) 

which are monotonously increasing functions for < u>3,ii>4 < 1. We can now reexpress 
the (N 3 LO) TPS of Rj(b;c) (in powers of b) as (N 3 LO) TPS in powers of u> 3 or to 4 , by 
simply using the power expansions of Eq. (f41~|) . The advantage of using this form of TPS 
for resummations lies in the fact that the convergence radius for the power series of Rj is 
now r Wa = 1/2 in the u> 3 -plane, and r W4 « 0.42 in the u> 4 -plane, and the circle of convergence 
reaches thus the first IR renormalon singularity &(«;&) = 2 (fc = 3,4), in contrast to the case 
of nontransformed b. In this way, the hindering influence of the UV singularities (negative 
6's) has been significantly weakened. The mapping u> 4 (6) apparently suppresses even more 
strongly the influence of the UV singularities than w 3 (b), but probably less strongly the 
influence of the NLO IR renormalon singularity (6 = 3). The final formula for r T in this 
formulation follows directly from the form 




1 , dbiw) 

dx — : — exp 
dw 



(p 1 (ml;a )+c) b{w) 



0o 



sin (nb(w)) 



b(w) (1-&H) (1-&H/2) (1-6H/3) (1-&H/4) 

) , (42) 



{cMw) /f3 ) Clb(w)/l30 e- cl6 H/A> _ 

x r(i+ Cl 6H//3 ) m*m(fi(w);£ 



where /•,(&) 's are fiven by Eq. (0), w stands for w 3 or w 4 , and we integrate actually only 
up to a u> max corresponding to b ~ 4 (for w 3 , 3 = 0.50536, w 3max = exp(— i0 3 ); for u> 4 , 
04 = 0.1 ~ 0, w 4m ax = exp(— i0 4 ) w 1). This upper bound on 6 is justified because the 
contributions from higher fe's are very strongly suppressed due to the exponent in Eq. ([12|). 
Strictly speaking, the path in the w 3 -plane should be along the positive axis (below it) from 
w 3 = to u> 3 = 1, and then along the arc (inner side) of the unit circle between w 3 = 1 and 
u; 3 = Wz(b = 4) = exp(— i</> 3 ), as shown in Fig. [IJ. However, for practical calculations, it is 
much more convenient to use the integration along the ray u> 3 = x exp(— i</> 3 ) (0 < x < 1), as 
denoted in Eq. (f42[) and shown in Fig. |l[ Both paths give the same answer, since the closed 
contour in Fig. [I] does not contain any singularities of the integrand. In the case of w 4 , the 
trick is the same, and we choose 04 = 0.1 instead of 04 = for the ray (see Fig. 0), in order to 
avoid any possible problems with numerical instability that would otherwise arise from the 
too extreme vicinity of the integration path to the possible singularities of the integrand. 

At this point, we can thus already regard the N 3 LO TPS of the function Rj(b(wk); c) = 
exp[c b(wk)/ Po]Rj(b(wk)', 0) (fc = 3,4) asaformof resummation of Rj, solely via the mappings 
([4"0|). Consequently, expressions ( f4"2"D evaluated using the aforementioned N 3 LO TPS of 
Rj(b(wk)',c) can be regarded as our resummed predictions of r T = r r (a ), being functions 
of the mentioned strong QCD couplant a . There is, however, an additional freedom of 
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choosing the value of the constant c in Eq. Since the available series of Rj(b(iVk)', c) is 

truncated, the results of Eq. ( |4"2"[ ) will have some unphysical dependence on the value of c. In 
one of our approaches, we will choose the latter value by the principle of minimal sensitivity 
(PMS), i.e., by the condition 

^#^ = 0, (43) 
oc 

when using for Rj(b(wk)) their N 3 LO TPS forms. 

We can, however, proceed in a different way. The N 3 LO TPS of Rj(b(wk);c) can be 
further resummed using Pade approximants (PA's) The authors of Ref. |59|] presented 

compelling arguments that combining the conformal transformations with the PA-type of 
resummations in general leads to significantly improved results, at least when a sufficient 
number of terms in the power expansion are known. Especially the diagonal or almost 
diagonal PA's [^,fj{|, in our case [2/1]^ (wp.) and [l/2]-^_. (u> fc ), may represent an efficient 
way of extending the applicability of expressions for Rj into the region sufficiently far away 
from the origin (quasianalytic continuation). 

However, there is a possibility that PA's do not lead to an improvement. This is some- 
times the case when the TPS in question is known to a relatively low order, e.g. to the N 3 LO. 
Since our available TPS's of Rj(wk, c) (A; = 3, 4) are known to the N 3 LO (provided a specific 
value of dij is taken), we have to find criteria for keeping or rejecting the resulting PA's. The 
PA [2/l]ft.(wk) c) predicts one real singularity of Rj(b(wk);c), and [1/2]^ [wk\ c) two singu- 
larities which can be either real or mutually complex conjugate. Physically, Rj(b(wk); c) has 
the strongest singularities (UV and IR renormalon singularities) at b(wk) — — 1, —2, . . . and 



b(wk) = 3,4, . . .. This means that in the case of conformal transformations (ID]), the singu- 
larities of PA's should preferably be at Wk values corresponding to b= —1 or 3: u>3° le = — 1 
or +1; il>4° 1c = — 1 or +0.6. We will include in our analysis the predictions with those PA's 
[2/1]^ (wk) whose poles satisfy the latter conditions. For that, we will use our freedom to 
adjust the value of the constant c in Eq. (f4"2|). On the other hand, PA's [1/2]^. (wf.; c), which 
have two poles, only rarely satisfy approximately the aforementioned conditions simultane- 
ously. 



Furthermore, in practical calculations, we prefer to use for the (|39|) in the Steven- 
son's equation fl3"8p the PA [2/3]m^g(a;), particularly because of its reasonable singularity 
structure (x po \ e ~ 0.311, corresponding to a^ ole ~ 0.98). The latter signals the breakdown of 
perturbative QCD (pQCD), and this choice of /3^(x) has been used previously by some of 
us in Refs. |p2] , pl|| . Since the effective energy in the discussed QCD observable r T is relatively 
low E ~ m T < 2 GeV, we can expect that this choice is not entirely irrelevant numerically. 



5 Variation of c in Eq. ( f42| ) corresponds to changing somewhat the numerical procedure used. 

6 PA [n/m]ji(w) to a function R(w) is a ratio of polynomials in w of degree n (numerator) and 
m (denominator). The power expansion of [n/m]n(w) must reproduce the terms of the power 
expansion of R(w) up to, and including, the term ~ty n+m . PA [n/m]n(w) can be determined if 
we know the TPS of R(w) up to, and including, the term ~ui n+m . PA's [n/n] are called diagonal. 
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Later we will comment on how much the results change when employing the (N 3 LO) TPS 
for /%jg(x), instead. 



V. PREDICTIONS FOR a s ; THEORETICAL UNCERTAINTIES ESTIMATE 

As discussed in the previous section, for a given choice of a = a s (m^; MS)/7r and of 
constant c, we calculate the observable r T = r T (a ;c) (|42"|), using for Rj(b(wk); c) one of the 
variants (^)- (|30|) (j = 1, 2, 3), and one of the two conformal transformations ( [411) (^ = 3,4), 
and for the resummation of Rj(b(wk)) either the N 3 LO TPS or [2/1] PA. We thus use, at 
given c, altogether 3 x 2 x 2 = 12 ways of calculating r T as function of ao = a s (m^., MS)/7r. 
The value c is adjusted as described at the end of the previous section: (a) according to 
the PMS (|3|) when using N 3 LO TPS for Rj(b(wk))] (b) according to the pole requirements 
w P° le = -1 or +1 « olc = -1 or +0.6) when using PA [2/l)^.(w k ) for R^w*)). We 

repeat the analysis for three choices of the N 3 LO coefficient (|33|) of the Adler function: 
d^ = 25, 15, 35 In the central case of cfjj = 25, and for r T close to the central measured 
value r T = 0.1960 flUf) , we display the relevant numerical results in Table |[ Displayed are 
the predictions of various approximants for the input values a s (m^.; MS) = 0.325 and 0.326. 
The first 14 entries are predictions of Eq. (f4*2"|) when using PA's [2/1]^ (wk) (j = 1,2,3; 
A; = 3,4), where c was adjusted so that these PA's yield the aforementioned pole values of 
the leading UV or subleading IR renormalon pole. The next six entries are predictions when 
using N 3 LO TPS for ^(6(w fc )) in Eq. © with c adjusted so that the PMS principle © 



is satisfied (local maximum). The last entry contains the arithmetic average f T of all 20 
predictions, for the aforementioned two values of a s (m^;MS). From that entry, we deduce 
that the central measured value r T = 0.1960 (|IBD is achieved by this arithmetic average f T 
at a s (m^;MS) = 0.3254. The uncertainty of the prediction due to resummation method 
( "truncation" error) can be estimated by comparing the aforementioned prediction with the 
prediction which differs the most from it, i.e., with the prediction using the PA [2/1]-^. (w^) 
with j = 3, k = 4 and c= 1.34 (see the 14th entry of Table |). This prediction differs from 
the aforementioned one by |5a s | tr . ~ 0.0024. 

Repeating the very same calculations (at the same values of c) for the correspondingly 
higher and lower input values of a s (m^;MS), we obtain the corresponding predictions of 
a s (m^.; MS) for the upper and lower bounds of the measured values of r T (|l8|) by demanding 
that the arithmetic average f T be equal to those upper and lower bounds. This procedure 
then results in the prediction 

a fS( m 2) = 0.3254 ± 0.0060 exp . ± 0.0060 EW ± 0.0052 CKM ± 0.0024 tr . (4° } = 25) . (44) 

It is gratifying that in the case of 4 0) = 25 so many PA cases give physically acceptable 
pole structure, and that the predictions for a^ s (m^) of the aforementioned 20 different ap- 
proaches differ from each other only little - they differ by at most 0.0024 from the prediction 
of their total average. If confining ourselves to just one conformal transformation, the cen- 
tral value changes by only ±0.0001 (0.3255 for k = 3; 0.3253 for k = 4). If using only the 
six entries from the approach with the PMS, the central value is 0.3259. If using only the 
entries with w\ oXe ~ —1, the central value is 0.3258. If using only the entries with w^ ole ~ +1 
and wf 016 ~ +0.6, the central value changes to 0.3244. 
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Further, for most of the approximants of the first eight entries of Table | (except those 
with j = 1), when using the same values of c but using PA's [1/2] instead of [2/1] for Rj(wk), 
the predictions r T differ from those of [2/1] by no more than 0.0005, and the predictions 
of a s {m 2 ) also by no more than 0.0005. In these cases of good agreement, the poles wl° le 
of the PA's [1/2] do not lie deep inside the &(wfc)-intervals [— 1,+3]. In most of the cases 
when the disagreement is larger, at least one of the poles of the [1/2] falls deep inside these 
intervals, often even within the intervals [—0.5, +2]. This offers us an additional evidence 
that our requirement that the poles of [2/1] be either at b(wk) = —1 (leading UV renormalon 
pole) of b(wk) = +3 (next-to-leading IR renormalon pole) was reasonable. 

In order to obtain the uncertainties of the prediction due to the uncertainty =25 ±10 
around the central value 4 = 25, we repeat the procedure for the case of 4 = 15 and 
4 0) = 35. 

In the case 4°^ = 15, the situation is very similar to the afore-discussed case of 4° )= 25. 
We obtain six entries when we require w^ ole — — 1: three entries each when w^ ole = +1 and 
^poie _ _|_q q- s j x en t r ies when applying the PMS condition (H3p with N 3 LO TPS for Rj(b(wk)) 



- see Table [J, when two choices a s (mJ;MS) =0.331 and 0.332 are made. The arithmetic 



average r T of all these 18 entries is then equal to the central measured value ( |18[) r T = 0.1960 



when a s (m 2 ; MS) = 0.3299, i.e., the value higher by 0.0045 than the corresponding central 
prediction of Eq. (|44]) of the case 4 0) = 25. 



In the case 4 = 35, the numerical situation is less favorable. None of the approaches 
with PA's [2/1]^. (wk) have acceptable solutions under the requirement u>£° le = —1. Such 
poles occur at c < —2, but the corresponding predictions of r T there are quite unstable under 
the variation of c. In two cases, acceptable solutions are obtained with the approach with 
PA's [2/1]^. (wk) when we require w^ ole = +1 or u>4° le = +0.6 - see the first two entries of 
Table fTJ. The approach with the PMS (^), when using N 3 LO TPS for Rj(b(wk)), appears 
to be more difficult as well; the derivatives dr T /dc are never zero, but are negative for all 
reasonable values of c; nonetheless, these slopes have the smallest negative values at specific 
values of c - see the corresponding six entries in Table |T|. Taking the arithmetic average 
f T of the eight entries of Table |TJ, we infer that the arithmetic average achieves the central 
measured value ( |T8| ) r T = 0.1960 for the value a s (m 2 ; MS) =0.3194, which is lower by 0.0060 
than the corresponding central prediction of Eq. fl44|) of the case 4 0) = 25. 



From the above considerations we infer that the uncertainty of the prediction of 
a s (m 2 ; MS) due to the uncertainty (|33"D of 4 is ±0.0060. 

There is yet another theoretical uncertainty involved in the prediction (|44]) , connected 
with the choice of the renormalization scheme (RSch). As shown in the Appendix, when 
taking a RSch different from the modified 't Hooft (mtH) RSch = c\ (k = 2, 3, . . .), the 
integral transformations involved become more complicated. Then we end up with, instead 
of the simpler formula (f42"D, the more complicated one ( |A.22j ). We note that the leading RSch 
parameter C2 has the values C2 ~ 3.16; 4.47; 6.58; 5.24 in the mtH, MS, principle of minimal 
sensitivity (PMS) @, and the effective charge (ECH) [[| RSch's (where we take ^/ = 3; and 
the PMS and ECH RSch's refer to the N 3 LO TPS of the Adler function). Therefore, we 
estimate as a characteristic deviation from C2 = c\ ~ 3.16 the value C2 ~ c\ + 3.4. We then 
perform the analysis in the RSch with the latter value of c<i (all the other c^'s unchanged), 
using formula ( |A.22| ). We use the RScl parameter value t; 2 =fi 2 /Q 2 = 1 [Q 2 = m 2 exp(iy)], 
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for df^ we use the central value c?!j = 25, and for the /3-functions we use again the [2/3] PA 
form. The results of the analysis are written in Table [TV]. The variation of the parameter c 
allowed us to obtain the desired locations of the poles only in a few cases. The arithmetic 
average values r T of Table |V| reach the central value prediction r T = 0.1960 when the 
the coupling is a s (m^; MS) = 0.3285, which is by 0.0031 higher than the central value in 
Eq. (44). Concerning the change of the next-to-leading RSch parameter C3, we note that 



c 3 « 5.6; 21.0; 36.8; 16.1 for the mtH, MS, PMS and ECH RSch's. Therefore, we choose a 
characterictic deviation of C3 from the mtH value C3 = cf ~ 5.6 to be C3 = cf+31. Completely 
analogous analysis as in the case of the C2-deviation leads to the results of Table |V[ The 
arithmetic average values f T in Table [V] reach the central value r T = 0.1960 for the coupling 
a s (m^.; MS) = 0.3274, which is by 0.0020 higher than the central value in Eq. ([0]). Therefore, 
adding the deviations 0.0031 and 0.0020 in quadrature gives the estimated uncertainty due 
to the changes of the RSch to be ±0.0037. 

We can thus add to the prediction (03) the discussed uncertainties due to the variation 



of and of the RSch, resulting in our final prediction 

af^iml) = 0.3254 ± 0.0060 cxp . ± 0.0060 EW ± 0.0052 CKM 

±0.0060 W3 ± 0.0037 RSc h ± 0.0024 t , , (45) 
= 0.3254 ± 0.0060 cxp . ± 0.0079 ew +ckm ± 0.0074 th . . (46) 

In the last line, we added the corresponding uncertainties in quadrature; the combined un- 
certainty due to the uncertainty of djj (|33|) , the resummation ("truncation") uncertainty, 
and the RSch uncertainty we call the theoretical (th.) uncertainty. This combined uncer- 
tainty is comparable with the two other uncertainties in Eq. (f46"|). If we use for the MS 
/3-function in Eq. (^) the N 3 LO TPS form instead of the [2/3]oj^g(x) PA form, we obtain, 
in a way completely analogous to that described in Table [| the central value prediction 
a s {m^ MS) =0.3245. This is by 0.0009 lower than the central value prediction in Eqs. (^5|)- 
([46|), indicating that those nonperturbative effects which originate from the behavior of the 
/3-function are not strong. The fact that we used in the Borel integral a finite b(w) mayL (=4) 
does not influence the results. Namely, if we increase this quantity to &(w) max = 5 [corre- 
sponding for w 3 to 03 = 0.64350 in Eq. fl42|), and for w 4 to 4 = 0.40272], the predictions for 
a s {m1) change by ~10~ 6 , i.e., insignificantly. 

We then RGE-evolved the result (UMll) f rom tne scale AJ = m r « 1.777 GeV to the 
scale M z = 91.19 GeV. We used again the aforementioned [2/3] ^^(x) PA form of the /%g- 



function, which is based on the known four-loop N LO TPS form of /3j^g |56]]. Therefore, we 



employed the corresponding three-loop matching conditions |)0[ for the flavor thresholds. 
The matching was performed at /i(iV/) = K,m q (Nf) with the choice k = 2, where n{Nf) is 
the scale above which Nf flavors are assumed active, and m q (Nf) is the running quark mass 
m q (m q ) of the iVy'th flavor. We further assumed m c (m c ) = 1.25 GeV and m^m^) = 4.25 
GeV |3[. We thus obtain from Eqs. (H)-© 

a MS(M|) = 0.1192 ± 0.0007 cxp . ± 0.0010 ew +ckm ± 0.0009 th . ± 0.0003 evol . , (47) 
= 0.1192 ±0.0015 . (48) 



In Eq. ([48]), we added all the uncertainties in quadrature. In Eq. fl47j) , we included the 
uncertainties due to the RGE evolution, which come primarily from varying n from 1.5 to 



17 



3, and from varying the quark masses m c (m c ) = 1.25 ± 0.10 GeV and m^m;,) =4.25 ± 0.15 
GeV (see Ref. [22] for more details). 

If we repeat the calculation of a s (m*) and a s (M 2 ) with the N 3 LO TPS function [in 
Eq. (|38| ) and in the RGE evolution from m 2 to M 2 ], instead of the used PA [2/3]gj^g(x), the 
central value prediction a s (M 2 ; MS) =0.1192 remains unchanged up to the displayed digits. 
This is so because this change of /3j^g causes the central value prediction of a s (m 2 ; MS) to be 
by 0.0009 lower [as already mentioned after Eqs. (|45| ) - (|46|) 1, but then the RGE evolution to 
/z 2 = M 2 with the changed fl^g pushes the result up, approximately neutralizing the former 
effect. 

In the analysis leading to the results P5j)-(j4"gp we assumed that the power-suppressed 
terms, apart from those from the quark masses, do not contribute to the considered observ- 
able R T , as already emphasized in Sec. [E|. As mentioned in that section, the inclusive (V+A) 
fit of the ALEPH Collaboration [p7], within their framework, predicted the contributions of 
the (massless) power-suppressed terms to the canonical observable r T to be consistent with 
zero: <5r Tj pg = 0.000 ± 0.004. If we assumed that the latter estimates were valid also in our 
framework,[] this 5r r ps would have to be subtracted from the values given on the right-hand 
side of flT8]), resulting in an additional, "PS" -uncertainty, term ±0.0040 for the r T . This 
would in turn give an additional approximate uncertainty ±0.0041pg in the result (|45|)-([46D 
for a s (m 2 ), and ±0.0005ps in the result ( $T\) for a s (M 2 ). The combined uncertainty ±0.0015 
for a s (M 2 ) in (|48|) would increase to ±0.0016. All the central value predictions in (|45|)-(f48l) 
would remain unchanged. 

There are at least two indications that the above results p5|)-(^6|) and (|47|)-(|4*8|) are 



not wrong. In the approach of Ref. [Effil , which involved ordinary Borel transform of D{Q ) 



and where the RScl was fixed according to the (local) PMS and we used the MS RSch, the 
resulting predictions were very similar: a s (m 2 ) = 0.3267 ± 0.0062 cxp ± 0.0082 ew+CK m ± 

0. 0073 mcth . and a s (M 2 ) = 0.1193 ± 0.0007 exp . ± 0.0010 ew +ckm ± 0.0009 meth . ± 0.0003 cvo i. 
= 0.1193 ± 0.0015.0 Especially the latter values are in virtual agreement with Eq. (|48|). 

Yet another indication that the results presented here are correct comes from repeating 
the entire resummations, but this time without employing the conformal transformations 
b = b(w). The /^-functions were again taken in the [2/3] PA form. We carried out the 
resummation with the Borel integration in Eq. (^) again up to 6 max = 4. For simplicity, we 
fixed this time the c value to c = 0. In each of the cases djj =25, 15, 35, we excluded from the 
analysis the approaches with those PA's which give physically unacceptable pole structure, 

1. e., which have poles well inside the 6-interval [—1,3]. Those of the PA's [2/1]^ (6) and 
[1/2]^ (6) which were not excluded do not contain poles in the 6-interval [—0.8,3]. 

In the case df^ = 25, the nonexcluded approaches were those with [2/1]^ (b) for j — 1, 2, 3; 
and [1/2]^. (b) for j = 1,2. Gratifyingly, the a s -predictions of these five approaches differ 



7 See the discussion in Sec. [0| about the differences between our and ALEPH's framework. 

8 The method (meth.) uncertainty in [^] is the combination of uncertainties from Sd^, the 
truncation (resummation), and the RScl and RSch ambiguities; the method uncertainty there thus 
corresponds to our theoretical (th.) uncertainty in Eqs. (|46|) and (0). 
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from each other only little: \/\a s {m? T )\ < 0.0004. We took the arithmetic average r T of all 
five r r -predictions, and obtained a s (ml) = 0.3257 ±0.0061 cxp ± 0.0061ew ± 0.0053ckm (for 
4 0) = 25). 

In the case d% = 35, the nonexcluded PA approaches were those with [2/1]^. (b) for 

j = 1,2. In addition, we did not exclude the approaches with the N 3 LO TPS of Rjip) for 
7 = 2,3, because their predictions are close to the aforementioned PA approaches. We then 
took the arithmetic average f T of all four r T predictions and obtained for the central value 
a s (m^) = 0.3191 [corresponding to the central value r T = 0.1960 in Eq. (0)], which is by 
0.0066 lower than the aforementioned central value 0.3257 of the d^ =25 case. 

In the case G?j$ = 15, the nonexcluded approaches were those with [2/1]^ (6) and 
[1/2]^ (6). Using the arithmetic average f T of these two r r predictions, we obtain the central 

value a s (m^) = 0.3293, which is by only 0.0036 higher than in the d^ = 25 case. 

This leads us to the following predictions of our method, when no conformal transfor- 
mation is used 

af^(m 2 T ) = 0.3257 ± 0.0061 cxp . ± 0.0081 ew +ckm ± 0.0066^ , (49) 
a MS(M z 2 ) = 0.1192 ± 0.0007 exp . ± 0.0010 ew +ckm ± 0.0008 5d3 ± 0.0003 evol . , (50) 



which is in almost complete agreement with the predictions (45) and (f|7]), obtained by 
the employment of the two conformal transformations ([4"0"D. In Eqs. ( f4"9"D and OHO]), we did 
not include the uncertainties due to the resummation (truncation) and due to the RSch 
ambiguities, because we regard these two predictions only as an additional cross-check of 
our main predictions and uncertainty estimates (f45|)-(f43). 



VI. SUMMARY 

We calculated the hadronic tau decay width r T by employing the contour integral form 



for this quantity in the complex momentum plane (16) and the modified Borel transform for 
the associated perturbative massless Adler function. By choosing a special renormalization 
scheme (modified 't Hooft scheme: Ck = c\, k = 2, 3, . . .), the integrand of the modified Borel 
transform is renormalization scale invariant. In our approach, we explicitly account for the 
structure of the leading infrared (IR) renormalon of the Adler function via the corresponding 
ansatz. Further, to accelerate the convergence, i.e., to minimize the resummation (trunca- 
tion) uncertainties, we employ two different conformal transformations which "map away" 
all the renormalon singularities, except the leading and subleading IR renormalons, onto 
the unit circle. The correct location of the leading ultraviolet (UV) renormalon or of the 
subleading IR renormalon is enforced by employing Pade approximants for the truncated 
perturbation series of the functions associated with the modified Borel transform. The Borel 
integration, in this appraoch, turns out to have suppressed renormalon ambiguity for r T at 
the leading IR renormalon singularity, and the ambiguity due to the subleading renormalons 
is strongly suppressed by the exponent in the Borel integral. We neglect in the observable 
r T all the possible power correction terms (except the dimension d= 4 quark mass terms), 
because the results of the ALEPH analysis [^7] suggest that such terms are consistent with 
zero or negligibly small even in our resummation framework. 
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Our analysis predicts the values of a^ s (m 2 ) and a^ s (M^) given in Eqs. d45|)-(f48D- These 
predictions agree well with the results obtained in our previous analysis |2lJ of r T where 
we employed the ordinary Borel transforms. The latter transforms have significantly dif- 
ferent behavior, expansions and the strengths of the renormalon singularities, and their 
integrands in the Borel integral are, in contrast to the present approach, renormalization 
scale dependent. Therefore, our present predictions represent a powerful reconfirmation of 



the predictions of |p2| . We consider this to be important, because analyses of r T which 
do not involve Borel transforms and do not account for the leading renormalon structure 
of the associated Adler functions ||2T|j6l] , [lO|Jl9| , |62l| give predictions for a s {m^) and a s (M^) 
which significantly differ from our predictions and significantly differ among themselves, as 



already emphasized in ||22|| . On the other hand, if accounting for the renormalon struc- 
ture via a large-/?o resummation of the ordinary Borel transform and employing an ECH- 
related resummation of r r , as performed by the authors of Ref. ||, their predicted values 



a 



MS (M Z 2 ) = 0.120 ± 0.002 come significantly closer to our prediction (|7|)-([4|). Further, 
our prediction a^ s (M^) = 0.1192 ± 0.0015 is completely compatible with the world average 



0.1184 ± 0.0031 as given in Ref. but somewhat less compatible with the world average 



0.1173 ± 0.0020 as given in Ref. g>3|. 
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Appendix A. MODIFIED BOREL TRANSFORMS IN THE GENERAL 
RENORMALIZATION SCHEME 

For a function f(y) with the (asymptotically divergent) expansion around y = 

oo 

f(y) = i + Y,fny n , (A.l) 

n=l 



the modified Borel transform J r f{Q was introduced by the authors of via the following 
expansion: 

OO -| 

TAz) = 1 + V f n7 r- r -z n , (A.2) 

A ; ^ {n + Cl z){n-l + Cl z) ■ ■ ■ (l+ciz) V ; 

where c\ is the coefficient at the two-loop term of the /^-function (JTTJ) . The authors of [12" 
further showed that the following integral transformation connects f(y) and J~f{z): 

f(y) = -(1 - cm) f dz e-/» ( V ~y iZ —±—F{z) . (A.3) 
y Jo \zj Y{l + ciz) 

Further, they showed that there corresponds to each singularity ~ (R — z)~ a ~ ClR of the 
ordinary Borel transform Ff(z) = 1 + J2fnZ n /n\ a singularity ~ (R — z)~ a of the modified 
Borel tranform J 7 f(z) 
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Ff{z) 



oo e 

(R - z)~ a {1 + O [(R-z) ln(l -*/#)]} 



(A.4) 



On the other hand, the perturbative expansion of the massless Adler function D(Q 2 ) has 
the canonical form (TO), with a = a(/z 2 ; c 2 , c 3 , . . .) being the QCD couplant, and the coeffi- 
cients dk having a specific RScl and RSch dependence dk = rffc(£ 2 ; c 2 , . . . , c*;) (£ 2 = fi 2 /Q 2 ) 
determined by the requirement of the RScl and RSch independence of D{Q 2 ) (cf. also 0, 
first entry) 



dt = 4 0) + A)ln£ 2 , 

d 2 = 4 0) + (di-d? 2 ) + ci(rfi-4 0) ) 

d 3 = 4 0) + ? > {d 1 d 2 -df ) df ) ) - 2(d\-df )z ) - ( Cl /2)(d\-df )2 ) 



( C2 -4 0) : 



+(c 2 d 1 -4 0) 4 0) )-(i/2)(c 3 -c: 



(0)v 



(A.5) 
(A.6) 

(A.7) 



d 4 = df + 4(d lC / 3 -< ; 4 u - ( Cl /3)(d 3 -4 u; ) + (s/3)(4-4 u;z ) 

-(28/3)(4^ 2 - 4° )2 4 0) ) + (2c 1 /3)(rf 1 rf 2 -4 0) 4° ) ) + (l/3)(c 2 rf 2 -4 0) rf : 



7(0)7(0) 



(o)^ 



j(0)2 N 



(Oh 



+ (14/3)(4-4° )4 ) - (4/3)(c 2 4-4 0) 4 0)2 ) + (c 3 ^i-4 0) 4 0) ) - (l/3)(c 4 -4 0) ) . (A.8) 



Here, the coefficients d^ are at the RScl jj, 2 = Q 2 (£ 2 = 1) and in the RSch c 
Comparing expansions (|I0|) and (|A.1| ), and inspecting the integral transformation ( |A.3| ), we 
may identify 



(o) jo) 

2 5 C 3 5 



a = y; D{Q 2 ) 



yf{y) 



(A.9) 



which immediately leads to the relations fi = d\—Ci, fk = d k — cidk-i {k > 2). If we use 
b = f3oz as the Borel variable, this allows us to write expansion (|A.2|) as 



^(6;e 2 ;c 2 ,c 3 ,...) = 1 + 



{dt-cx) 



+ 5n 



(l+ Cl 6//3 ) \Po, 
(d n — Cid n _i) 



n=2 



+ c 1 b/(3 )(2+c 1 b/(3 ) ■ ■ ■ (n + db/Po) 



and the integral transformation ( |A.3|) as 



1 POO 

D(Q 2 ) = — / dbexp 

1 f°° 
— — \ dbexp 
Po Jo 



-ab/Po 



1 



-— - + ci ln(cia; 



r(i+ Cl 6/A)) 

kTo) 



(A.10) 



r(i+ Cl 6/A) 



-T D {b) 



(A.ll) 
(A.12) 



If the ordinary Borel transform D(b; £ 2 ; c 2 , . . .) (|2T|) has a singularity of the form ~ — 
b/(3o)~ a ~ ClR , then ~ ( J R-z)- Q ~ Cli? [l + C(( J R-^) 1 )], due to the simple relation dF f /dz = 
dD/dz — C\D. Therefore, we can write the singularity relation between the ordinary Borel 
transform D and the modified Borel transform To^Tj in complete analogy with Eq. ( |A.4|) 
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D(b;e-,C2,...)~(K-b)- a - ClK/l30 
T D {b- i 2 - 02, . . .) ~ {K-b)- a {1 + O [(«-&) 1ii(1-&/k)]} 



(A.13) 



where k = /5o-R. We may use the subtracted Stevenson equation [cf. Eq. (|38])] to reex- 
press (partly) the expression (l/a + c\ ln(cia)) in the exponential of Eq. ( |A.12| ) in terms of 
ln(/i 2 /Q 2 )=ln(^ 2 ) and of the invariant pi(Q 2 ) (H) 



- + Cl ln(cia) = pi (Q 2 ) + 4 0) + A) In f - f 
a jo 



.IT 



(3(x) 



(A. 14) 



where /3(x) is in the RSch considered, i.e., its expansion around x = is: f3(x\ c 2 , . . .)//3 : 
— a; 2 (l + Cia;+C2a; 2 + - • ■). Therefore, Eq. ([A. 12 ) can be rewritten 



D(Q 2 )=D a(W;c 2 ,...);e 2 ;c 2 



1 r 00 

— / <i&exp 
Pn Jo 



(pi(g 2 )+/3olna& 

A) 



x exp 



b r a 



Pa Jo 



dx 



(1 — CiX) 



x- 



+ 



( Cl 6//5 ) Cl6//3 ° 

x r(i+ Cl 6//3 ) exp 



4^6 

' A) 



/5(x;c 2 ,...) / 
FD(b;£ 2 ;c 2 , 



(A.15) 



In the general RSch (c 2 , c 3 , . . .) and at the general RScl £ 2 = fi 2 /Q 2 , the modifed Borel 
transform function Tr>{b\ £ 2 ; c 2 , . . .) is related to the function J-£>(b; 1; c 2 , . . .) in the following 
way: 



e- 6ta ^(6; £ 2 ; c 2 , c 3 , . . .) = T D {b; 1; cf, of, . . .) - i(c 2 -c 2 0) ) 



1 

+ 12 



+ 72L 



2(c 2 -c 2 )(/? lne 2 ) + (c 2 -4 0) )(ll Cl -44 u O - (c 3 -4 u 



j(0h 



,(0h 



3(c 2 -^)(/5 ln^)(24 0) -/5 lne 2 ) + (22c? + 3c 3 -25 Cl c 2 )(A> lnf 



-(c 2 -4 0) )(85c 2 -4c 2 + 5c 2 u; -50 Cl < ; + 94 Uj ) + (c 3 -c 3 Uj )(13 Cl -3dn - (c 4 -c. 
+0(6 5 ) . 



m , O j(o) 



(oh 



j(0h 



,(0h 



(A.16) 



This can be shown, for example, by using transformation formulas (|A.5| )- (|A.8| ) in the expan- 
sion ( |A.10|) . The above formulas (|A.16|) and ( |A.15 ) show that the modified 't Hooft (mtH) 
scheme (ck = c\, k — 2, 3, . . .) is a remarkable RSch choice in the discussed Borel transforms 



e- feln ^(6; e; 4, 4, . . .) = F D (b; 1; c 2 , 4, . . .) 



[1-ctx) 



+ 



00 



x- 



(3(x; 4,4,...) 
and thus Eq. ( |A.15| ) reduces in this case to 



(A.17) 
(A.18) 
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D(Q 2 ) = D(a(eQ 2 ;clcl..y,e;cicl 



Po 




(ci&/A)) cl6/A) 

x r(i+ Cl 6//5 ) exp 



This is just the integral transformation fl23|) , with expansion 
(3 \n^ 2 [note that d x = df } = di(£ 2 = 1)] 



Po 

...). (A.19) 
)f) and the constant c = 



D(6;c) =e a//3o D(6;0) 



(ci6//3b 



\cib/po 



T(l + Cl b/p 



■ exp 



(c -4 0) )& 

Po 



F D (b;l;4,4,...) 



(A.20) 



The above expression ( |A.19| ) shows the remarkable property of the mtH RSch: the whole 
integrand in the modified Borel transformation in the mtH RSch is RScl-independent (£ 2 - 
independent). This appears the main reason why Grunberg [|T called this Borel transforma- 
tion "renormalization scheme invariant." Strictly speaking, we see that this transformation 
may be called "renormalization scale (RScl) invariant," with the choice of the mtH renor- 
malization scheme (RSch). On the other hand, if we abandon the mtH RSch in the general 
class of Borel transformations ( |A.15| ), the integrand becomes explicitly RScl-dependent 



d [integrand (b; £ 2 ; c 2 , . . .)] 
P^l 2 

oc ( Cl 6//3 ) Cl6//3 ° I I - I 1 
= ( Cl 6//5 ) Clb//3 ° 



b\ 




Po) 




b\ 




Po) 


\ 



{l-c x a)P{a] c 2 ,...) 



+ 0{b 2 



p a 2 

1 — (1 — Cia)(l+Cia+c 2 a 2 -|- ■ ■ 



0{b 2 ) 



(A.21) 



In order to obtain the expression for r T in terms of the modified Borel transform 
To{b'i £ 2 ; c 2 , . . .) in the general RSch, the contour integration fliTf ) in the complex momentum 
plane has to be performed on the massless Adler function expression ( |A.15| ). Further, we 
can perform in addition the conformal transformation b = b(w) of the types (|4lD and the ray 
integration trick in the u>-plane as explained in Figs. 0,0. The procedure is analogous to 
the procedure leading to formula (42) in the mtH RSch, and we end up with the following 
formula in the general RSch: 



27r/3 



-Re<^ e" i9 



o 



1 , db(w) 
ax — - — exp 

aw 



(pi(m 2 ; a ) + c) b(w) 



x 



( Cl b(w)/p ) clbM/l30 e~ c ^l^ 
T(l + Cl b(w)/p ) 



Po 



x exp 



1-&H/2) 

b(w) / ,a (? 2,7i ? ex p( i 2/); c 2r-) 



dy{l + e iy f{l-e iy )e- ih{w)y 



Po Jo 



dx 



[l-cix) 



X* 



+ 



Po 



P(x;c 2 , 



(A.22) 



where the weight functions fj are given in Eq. (|37D and the Borel transform functions Rj 
are defined in analogy with Eqs. (^7|)-(pC|) 
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i2i(6;c;e 2 ;c2,...) = (1-6/2) exp[(£+ Cl -4 0) )6//3 ]e- 61n «V D (6; c 2 , c 3 , .. .) , (A.23) 
R 2 (b; g; e 2 ; c 2 , . . .) = Si(6; c; £ 2 ; c 2 , . . /n j^ffi > (A.24) 

75 3 (fc; c; ft c 2 , . . .) = fl l( 6; c; ft c 2 , . . ) {l + f^f+f M . (A.25) 

(3 + c^/Po) 2 

We introduced an additional freedom factor exp(— c 6//?o) in the exponential in Eq. ( |A.22[) , 
which is then offset by the factor exp(+c b//3o) in the functions Rj, in analogy with the case 
of the mtH RSch Eq. fl42|). The functions Rj(b(w); £ 2 ; c; c 2 , . . .) can be resummed, either as 
PA [2/1]^ (w) or as simple N 3 LO TPS. The function exp(-6 In ^)T D {b] £ 2 ; c 2 , c 3 , . . .), which 

is the source of the RScl-dependence in R^s, has only a weak RScl-dependence, as can be 
seen from Eq. (|A.16|) 

e- hh ^F D (b- e; c 2 , c 3 , . . .) = T D {b- 1; c 2 , c 3 , . . .) + 

\^- c ^{j}j (/W 2 ) + C>(& 4 ) • (A.26) 

We further see from Eq. ( A.22|) that in the general RSch we cannot perform the contour 
integration over dy analytically, in contrast to the mtH RSch Eq. 



24 



REFERENCES 



[1] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983). 

[2] P. M. Stevenson, Phys. Rev. D 23, 2916 (1981); Phys. Lett. 100B, 61 (1981); Nucl. 

Phys. B203, 472 (1982). 
[3] G. Grunberg, Phys. Lett. 95B, 70 (1980), HOB, 501(E) (1982); 114B, 271 (1982); 

Phys. Rev. D 29, 2315 (1984). 
[4] A. L. Kataev, N. V. Krasnikov and A. A. Pivovarov, Nucl. Phys. B 198, 508 (1982) 

[Erratum-ibid. B 490, 505 (1982)] ||hep-ph/9612"326| . 
[5] A. Dhar and V. Gupta, Phys. Rev. D 29, 2822 (1984); V. Gupta, D. V. Shirkov, and 

O. V. Tarasov, Int. J. Mod. Phys. A 6, 3381 (1991). 
[6] H. Lu and S. J. Brodsky, Phys. Rev. D 48, 3310 (1993) [|hep-ph/9211254|1 ; S. J. Brodsky 
and H. J. Lu, ibid. 51, 3652 (1995) |hep-ph/94052T8; l; S. J. Brodsky, G. T. Gabadadze, 
A. L. Kataev and H. J. Lu, Phys. Lett. B 372, 133 (1996) ||hep-ph/951236ll ; S. Groote, 
J. G. Korner, A. A. Pivovarov and K. Schilcher, Phys. Rev. Lett. 79, 2763 (1997) pep 



ph/9703208|| ; S. J. Brodsky, J. R. Pelaez and N. Toumbas, Phys. Rev. D 60, 037501 



(1999) ||hcp-ph/98lG4^4] . M. Binger, C. Ji, and D. G. Robertson, ibid. 61, 114011 (2000) 



||hep-ph/991224T 



[7] C. J. Maxwell, |iep-ph/9908"463l C. J. Maxwell and A. Mirjalili, Nucl. Phys. B577, 209 

(2000) ||hep-ph/ 0002201 . 
[8] C. J. Maxwell and D. G. Tonge, Nucl. Phys. B 535, 19 (1998) ||hep-ph/97053l4 |. 
[9] C. J. Maxwell and A. Mirjalili, |hep-ph/0103T64l . 
[10] J. G. Korner, F. Krajewski, and A. A. Pivovarov, Phys. Rev. D 63, 036001 (2001) 
||hep-ph/0002T66l . 

[11] D. V. Shirkov and I. L. Solovtsov, Phys. Rev. Lett. 79, 1209 (1997) [|hep-ph/9704333|| ; 
K. A. Milton, I. L. Solovtsov, and O. P. Solovtsova, Phys. Lett. B 439, 421 (1998) 
hep-ph/98095T0l ; Phys. Rev. D 64, 016005 (2001) ||hep-ph/0102"25l ; Eur. Phys. J. C 



[12 
[13 

[14 

[15 

[16 
[17 



14, 495 (2000) ||hep-ph/000303T| ; I. L. Solovtsov and D. V. Shirkov, Phys. Lett. B 442, 
344 (1998) ||hep-ph/971 12511 ; Theor. Math. Phys. 120, 1220 (1999) ||hep-ph/9909305 
D. V. Shirkov, ibid. 127, 409 (2001) ||hep-ph/0012"283T 



[18 

[19 
[20 



D. S. Kourashev, |hep-ph/99124TUt |hep-ph/0010072 . 

S. J. Brodsky, E. Gardi, G. Grunberg, and J. Rathsman, Phys. Rev. D 63, 094017 
(2001) ||hep-ph/ 00020651 . 

I. L. Solovtsov, Phys. Lett. B 327, 335 (1994); 340, 245 (1994); K. A. Milton, I. L. 
Solovtsov, and O. P. Solovtsova, Eur. Phys. J. C 13, 497 (2000) ||hep-ph/0003Tm| . 
T. Lee, Phys. Re v. D 56, 1091 (1997) |hep-th/96110T0| ; Phys. Lett. B 462, 1 (1999) 
||hep-ph/9908225l . 



I. Caprini and J. Fischer, Phys. Rev. D 60, 054014 (1999) ||hep-ph/9811367|] . 

M. A. Samuel, J. Ellis, and M. Karliner, Phys. Rev. Lett. 74, 4380 (1995) fhepj 

ph/950341l| ; J. Ellis, E. Gardi, M. Karliner, and M. A. Samuel, Phys. Lett. B 366, 



268 (1996) ||hep-ph/95093lH ; Phys. Rev. D 54, 6986 (1996) ||hep-ph/9607404| . 
E. Gardi, Phys. Rev. D 56, 68 (1997) ||hep-ph/96ll45l ; S. J. Brodsky, J. Ellis, E. Gardi, 
M. Karliner, and M. A. Samuel, ibid. 56, 6980 (1997) [ |iep-ph/9706467j . 
T. G. Steele and V. Elias, Mod. Phys. Lett. A 13, 3151 (1998) ||hep-ph/9902"2T7[ . 
G. Cvetic, Nucl. Phys. B 517, 506 (1998) ||hep-ph/9711406H ; Phys. Rev. D 57, R3209 



25 



(1998) [|hep-ph/9711487|| ; G. Cvetic and R. Kdgerler, Nucl. Phys. B 522, 396 (1998) fEep 



ph/9802248|| ; G. Cvetic, Nucl. Phys. B (Proc. Suppl.) 74, 333 (1999) [|hep-ph/9808273 



[21] G. Cvetic and R. Kdgerler, Phys. Rev. D 63, 056013 (2001) [|hep-ph/0006098|| ; G. Cvetic, 
Phys. Lett. B 486, 100 (2000) ||hep-ph/ 0003121 . 

G. Cvetic and T. Lee, Phys. Rev. D 64, 014030 (2001) ||hep-ph/010129T |. 
M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979); 
ibid 448 (1979); ibid 519 (1979). 

E. Braaten, S. Narison, and A. Pich, Nucl. Phys. B 373, 581 (1992). 
Y. L. Dokshitzer, G. Marchesini, and B. R. Webber, Nucl. Phys. B 469, 93 (1996) fhepj 



ph/9512336|| ; E. Gardi and G. Grunberg, JHEP 9911, 016 (1999) ||hep-ph/9908458 



G. Grunberg, ibid. 9811, 006 (1998) ||hep-ph/9807494 |. 
T. Lee, |hep-ph/0104306| . 

R. Barate et al. [ALEPH Collaboration], Eur. Phys. J. C 4, 409 (1998). 

R. Barate et al. [ALEPH Collaboration], Eur. Phys. J. C 11, 599 (1999) |Eep- 

ex/ 99030 15|| ; ibid. 10, 1 (1999) ||hep-ex/99030T4 |. 



M. Davier, [hep-ex/99120l4l . 

G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C 13, 197 (2000) fiepj 
ex/ 99080 13|| . 



S. J. Richichi et al. [CLEO Collaboration], Phys. Rev. D 60, 112002 (1999) |iepj 
ex/9810026|| . 



M. Davier, private communication; and talk given at the Sixth International Workshop 
on Tau Lepton Physics (Tau'2000 Workshop), September 2000, Victoria, Canada (to 
appear on [http:/ /tau2000.phys.uvic.ca/| ) . 

D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 
C. S. Lam and T. M. Yan, Phys. Rev. D 16, 703 (1977). 

K. Schilcher and M. D. Tran, Phys. Rev. D 29, 570 (1984). 
S. Narison and A. Pich, Phys. Lett. B 211, 183 (1988). 

E. Braaten, Phys. Rev. Lett. 60, 1606 (1988); Phys. Rev. D 39, 1458 (1989). 

A. Pich, Hadronic tau-decays and QCD, in Proc. Workshop on Tau Lepton Physics, 
Orsay, France, September 1990, eds. M. Davier and B. Jean-Marie (Ed. Frontieres, 
Paris, 1991), p.321 (see also preprint CERN-TH-5940-90, KEK Library). 
A. A. Pivovarov, Z. Phys. C 53, 461 (1992) [Yad. Fiz. 54, 1114 (1991), Sov. J. Nucl. 
Phys. 54, 676 (1992)]; F. Le Diberder and A. Pich, Phys. Lett. B 286, 147 (1992). 
A. Pich and J. Prades, JHEP9806, 013 (1998) ||hep-ph/ 9804352] ; ibid. 9910, 004 (1999) 
||hep-ph/ 9909244| . 



G. Grunberg, Phys. Lett. B 304, 183 (1993). 

L. S. Brown and L. G. Yaffe, Phys. Rev. D 45, 398 (1992); L. S. Brown, L. G. Yaffe 

and C. Zhai, ibid. 46, 4712 (1992) ||hep-ph/9205213|] . 

E. Braaten and C. Li, Phys. Rev. D 42, 3888 (1990). 

W. J. Marciano and A. Sirlin, Phys. Rev. Lett. 61, 1815 (1988). 

K. G. Chetyrkin, S. Narison, and V. I. Zakharov, Nucl. Phys. B 550, 353 (1999) fFiepj 
ph/ 98 112751 . 



S. J. Huber, M. Reuter, and M. G. Schmidt, Phys. Lett. B 462, 158 (1999) flnej> 



ph/9906358 



C. A. Dominguez and K. Schilcher, Phys. Rev. D 61, 114020 (2000) [ |hep-ph/9903385 



26 



K. Maltman, Phys. Lett. B 440, 367 (1998) ||hep-ph/9901239| . 
C. A. Dominguez and K. Schilcher, Phys. Lett. B 448, 93 (1999) ||hep-ph/98 11261 . 
K. G. Chetyrkin, A. L. Kataev, and F. V. Tkachev, Phys. Lett. B 85, 277 (1979); 
M. Dine and J. Sapirstein, Phys. Rev. Lett. 43, 668 (1979); W. Celmaster and R. J. Gon- 
salves, ibid. 44, 560 (1980). 

S. G. Gorishny, A. L. Kataev and S. A. Larin, Phys. Lett. B 259, 144 (1991); L. R. Sur- 
guladze and M. A. Samuel, Phys. Rev. Lett. 66, 560 (1991) [Erratum-ibid. 66, 2416 
(1991)]. 

M. Beneke and V. M. Braun, Phys. Lett. B 348, 513 (1995) ||hep-ph/ 94 112291] ; P- Ball, 
M. Beneke, and V. M. Braun, Nucl. Phys. B 452, 563 (1995) ||hep-ph/9502"300 |. 
M. Neubert, Nucl. Phys. B 463, 511 (1996) ||hep-ph/9509432|] . 
C. N. Lovett-Turner and C. J. Maxwell, Nucl. Phys. B 452, 188 (1995) 



ph/9505224 



A. H. Mueller, Nucl. Phys. B 250, 327 (1985). 

T. van Ritbergen, J. A. Vermaseren, and S. A. Larin, Phys. Lett. B 400, 379 (1997) 
||hep-ph/970T390l . 

G. Altarelli, P. Nason, and G. Ridolfi, Z. Phys. C 68, 257 (1995) ||hep-ph/950i240|1 . 
George A. Baker, Jr. and Peter Graves-Morris, Pade Approximants , 2nd edition, (En- 
cyclopedia of Mathematics and Its Applications, Vol. 59), edited by Gian-Carlo Rota 
(Cambridge University Press, 1996). 

U. D. Jentschura and G. Soff, J. Phys. AA 34, 1451 (2001) ||hep-ph/0006089| . 



K. G. Chetyrkin, B. A. Kniehl, and M. Steinhauser, Phys. Rev. Lett. 79, 2184 (1997) 
||hep-ph/9706430l . 

S. Groote, J. G. Kdrner, and A. A. Pivovarov, Mod. Phys. Lett. A 13, 637 (1998) 
||hep-ph/97032B8~ |. 

B. V. Geshkenbein, B. L. Ioffe, and K. N. Zyablyuk, Phys. Rev. D (to be published) 

||hep-ph/ 01040481 . 

S. Bethke, J. Phys. G G26, R27 (2000) ||hep-ex/000402TH . 

I. Hinchliffe and A. V. Manohar, Ann. Rev. Nucl. Part. Sci. 50, 643 (2000) flhel> 
ph/0004186|| . 



27 




w-plane 



w(b=3)=l 



FIG. 1. Integration in the W3-plane along the ray w 3 = xexp(— i(f> 3 ) (0 < x < 1, (f>s — 0.50536) gives 
the same result as the integration parallel to the positive real axis (0 < w < 1) and arc w — exp(— i(f>') 
(0 < <j>' < fo). 
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(j, k); approximant 


r T [a s (m 2 T ) = 


= 0.325; 0.326] 


c 


comments 


(1,3) 


[2/11 


0.19543 


0.19641 


+0.28 


w po\e _ _l 


(1,4) 


[2/1] 

L / J 


0.19524 


0.19621 


+0.17 


^polc ra _ L 


(2,3) 


[2/1] 


0.19533 


0.19631 


-0.27 


^polc w _ L 


(2,4) 


[2/1] 


0.19531 


0.19629 


-0.28 


^polc ra _J 


(3,3) 


[2/1] 


0.19570 


0.19668 


-0.15 


w Polc ~ _! 


(3,4) 


[2/1] 


0.19554 


0.19651 


-0.22 


^polo ra _ L 


(2,3) 


[2/11 

L / J 


0.19470 


0.19567 


-0.71 


^polc „ _ L 


(2,4) 


[2/1] 


0.19466 


0.19563 


-0.89 


^polc w _ L 


(1,3) 


[2/1] 


0.19503 


0.19600 


+1.95 


w Polo ~ +1 _ 


(1,4) 


[2/1] 


0.19587 


0.19686 


+2.275 


^polc ^ + Q g 


(2,3) 


[2/1] 


0.19648 


0.19748 


+1.23 


^polc _ +1 _ 


(2,4) 


[2/1] 


0.19746 


0.19847 


+1.61 


^pole _ + q g 


(3,3) 


[2/1] 


0.19686 


0.19786 


+0.97 


^pole „ +L 


(3,4) 


[2/1] 


0.19792 


0.19895 


+1.34 


^pole _ +Q g 


(1,3) 


TPS 


0.19449 


0.19546 


+0.92 


local max. 


(1,4) 


TPS 


0.19433 


0.19529 


+0.54 


local max. 


(2,3) 


TPS 


0.19545 


0.19643 


+0.47 


local max. 


(2,4) 


TPS 


0.19516 


0.19614 


+0.12 


local max. 


(3,3) 


TPS 


0.19568 


0.19666 


+0.31 


local max. 


(3,4) 


TPS 


0.19535 


0.19633 


-0.01 


local max. 


arithm. average r T 


0.19560 


0.19658 



TABLE I. The results r T of calculations according to Eq. (p2h , using 



transformations (|4l| ) (fc = 3,4) and the resummations of Rj(wk) (f 
(N 3 LO) TPS. The constant c was varied so as to achieve either 



w 



pole 



25, employing the conformal 
Q (j = 1, 2, 3) as PA [2/1] and as 



-1 or ai 



pole 



+1 K olc = +0.6) 

when using PA's [2/1], or stationarity according to the PMS condition (43) when using (N 3 LO) TPS. The 
results r T are given for the choices a s (m^;MS) = 0.325 and 0.326. The central measured value r T = 0.1960 
( |j~8| ) is achieved by the arithmetic average r T at a s {w? T \ MS)=0. 3254. 



29 



(j, k); approximant 


t \rv (m ^ 1 = 


= 331- 3321 


(I 


rom in pnts 


(1,3) 


[2/1] 


U.lyozo 


u.iy f i { 


i 1 on 


„,,P01C 1 


(1,4) 


[2/1] 


U.lyoUU 


u.iyoyz 


i i in 
+ 1.111 


„,,polc ^ 1 

w ~ — 1. 


(2,3) 


[2/1] 


U.ly ( Yl 


u.iyoUo 


i n qq 


„,,polc ^ 1 


(2,4) 


[2/1] 


u. ±yuoo 


u. iy i oi 


I U. i Q 


,„pole ^ I 


(3,3) 


[2/1] 


0.19776 


0.19871 


+0.55 


w Polc ~ _! 


(3,4) 


[2/1] 


0.19750 


0.19844 


+0.46 




(1,3) 


[2/1] 


0.19601 


0.19693 


+2.52 


w Polc w +1 _ 


(1,4) 


[2/1] 


0.19687 


0.19781 


+2.83 


^polc _ + Q g 


(2,3) 


[2/1] 


0.19758 


0.19853 


+1.98 


w Pole +1 _ 


(2,4) 


[2/1] 


0.19856 


0.19953 


+2.30 


^pole _ +0 _g 


(3,3) 


[2/1] 


0.19831 


0.19927 


+1.64 


^polc w +1 _ 


(3,4) 


[2/1] 


0.19942 


0.20040 


+1.97 


^polc _ + Q g 


(1,3) 


TPS 


0.19543 


0.19634 


+1.6 


local max. 


(1,4) 


TPS 


0.19525 


0.19616 


+1.3 


local max. 


(2,3) 


TPS 


0.19666 


0.19759 


+1.2 


local max. 


(2,4) 


TPS 


0.19640 


0.19733 


+0.9 


local max. 


(3,3) 


TPS 


0.19732 


0.19826 


+0.9 


local max. 


(3,4) 


TPS 


0.19703 


0.19797 


+0.6 


local max. 


arithm. average r T 


0.19702 


0.19796 



TABLE II. Results analogous to those of Table |. but for the case <4 = 15. The results r T are given 
for the choices a s (m^; MS) = 0.331 and 0.332. The central measured value r T — 0.1960 (|lj) is achieved by 
the arithmetic average f T at a s (rr? T \ MS) = 0.3299. 
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(j, k); approximant 


r T [a a (ml) = 


= 0.331; 0.332] 


c 


comments 




(1,3) 


[2/1] 


0.19716 


0.19823 


+1.02 


^pole 


Pa +1. 




(1,4) 


[2/1] 


0.19780 


0.19888 


+1.37 


^polc 


« +0.6 




(1,3) 


TPS 


0.19744 


0.19852 


+0.025 


dr T /dc pa 


-2.66 • 10' 


-3 


(1,4) 


TPS 


0.19645 


0.19751 


+0.0 


dr T /dc pa 


-2.92 • 10" 


-3 


(2,3) 


TPS 


0.19603 


0.19708 


+0.375 


dr T /dc pa 


-3.18 • 10" 


-3 


(2,4) 


TPS 


0.19558 


0.19662 


+0.075 


dr T /dc pa 


-3.09 • 10- 


-3 


(3,3) 


TPS 


0.19651 


0.19757 


-0.19 


dr T /dc pa 


-2.37-10- 


-3 


(3,4) 


TPS 


0.19602 


0.19707 


-0.49 


dr T /dc pa 


-2.29 • 10" 


-3 


arithm. average r T 


0.19662 


0.19769 



TABLE III. Results analogous to those of Table ffl, but for the case d 3 =35. The results r T are given 
for the choices a s (m 2 T ; MS) = 0.320 and 0.321. In the approach with the (N 3 LO) TPS of Rj(w k ), the PMS 
condition ( fllij ) is never exactly satisfied; there, the values of c were chosen so that the (negative) slope dr T /dc 
is the least steep. The central measured value r T = 0.1960 ( |l8|) is achieved by the arithmetic average f r at 
a s (m2;MS) = 0.3194. 



(j, A;); approximant 


r T [ot s (m 2 T ) = 


= 0.329; 0.330] 


c 


comments 


(2,4) 


[2/1] 


0.19521 


0.19613 


+0.01 


^pole w _ L 


(3,3) 


[2/1] 


0.19656 


0.19750 


+0.707 


^polc _ +1 


(3,3) 


[2/1] 


0.19770 


0.19866 


+1.25 


^pole „ +L 


(2,3) 


TPS 


0.19668 


0.19763 


+1.1 


local max. 


(2,4) 


TPS 


0.19633 


0.19727 


+0.8 


local max. 


(3,3) 


TPS 


0.19672 


0.19767 


+1.05 


local max. 


(3,4) 


TPS 


0.19638 


0.19732 


+0.75 


local max. 


arithm. average r T 


0.19651 


0.19745 



TABLE IV. Results analogous to those of Table [j], but for the different RSch: c-i = cf + 3.4, Cfe = 
(fc = 3, 4, . . .). Expression (A. 22) is used to calculate r T . RScl parameter is £ 2 = /j, 2 /Q 2 — 1 [Q 2 — m 2 exp(iy)]; 

=25. The results r T are given for the choices a s (m 2 ; MS) = 0.329 and 0.330. The central measured value 
r T = 0.1960 (|lq) is achieved by the arithmetic average f T at a s (m 2 ; MS) = 0.3285. 
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(j, k); approximant 


r T [a s (m 2 T ) = 


= 0.327; 0.328] 


c 


comments 




\9 /ll 
l z / L l 


0.19478 


0.19570 


+1.50 


^polc w _ L 


U, 4 J 


\0 /ll 


0.19453 


0.19545 


+1.40 






F2/11 


0.19582 


0.19676 


+1.11 


^polc w _ L 


(2,4) 


[2/1] 


0.19558 


0.19652 


+1.01 


w Polo w _ L 


(3,3) 


[2/1] 


0.19654 


0.19749 


+0.80 


^polc ra _J 


(3,4) 


[2/1] 


0.19628 


0.19723 


+0.70 


w Polo ra _ L 


(1,3) 


[2/1] 


0.19456 


0.19549 


+2.77 


w polc + L 


(1,4) 


[2/1] 


0.19538 


0.19632 


+3.07 


^polc pa +0 _g 


(2,3) 


[2/1] 


0.19612 


0.19708 


+2.25 


w po1c p. +1. 


(2,4) 


[2/1] 


0.19708 


0.19805 


+2.56 


^polc pa +0 6 


(3,3) 


[2/1] 


0.19695 


0.19792 


+1.90 


^polc pa +L 


(3,4) 


[2/1] 


0.19796 


0.19894 


+2.205 


^pole pa +Q g 


(1,3) 


TPS 


0.19400 


0.19492 


+1.9 


local max. 


(1,4) 


TPS 


0.19383 


0.19475 


+1.5 


local max. 


(2,3) 


TPS 


0.19528 


0.19622 


+1.45 


local max. 


(2,4) 


TPS 


0.19505 


0.19598 


+1.1 


local max. 


(3,3) 


TPS 


0.19602 


0.19696 


+1.1 


local max. 


(3,4) 


TPS 


0.19576 


0.19670 


+0.8 


local max. 


arithm. average r T 


0.19564 


0.19658 



TABLE V 

(A; = 2,4,5,...). Expression (A. 22) is used to calculate r 



Results analogous to those of Table IV, but for the RSch: C3 

RScl parameter is £ 2 



4 



1; dg" = 25 



results r T are given for the choices a s (m^;MS) = 0.327 and 0.328. The central measured value r T 
(O) is achieved by the arithmetic average f T at a s (m^; MS) = 0.3274. 



31.0, c fe = cf 

(0) _ OK Xhc 

0.1960 
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